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is solved, an adiabatic evolution is defined and the corresponding adiabatic condition is given. 
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I. INTRODUCTION 

The evolution of a system that interacts with its sur- 
rounding environment is fully given by a dynamical map 
that corresponds to a quantum stochastic process. The 
state, the environment and their correlations change 
with time. If the environment is assumed not to re- 
act on the system with memory, the Markov approx- 
imation can be taken in which these correlations are 
discarded to derive the Kossakowski-Lindblad master 
equation [l], 0, This theory extends quantum me- 
chanics beyond Hamiltonian dynamics, and find pow- 
erful applications in quantum opticsk I and quantum 
information Q . Many approaches [g, 0, IsM 03 HH EH 
have been proposed to solve this Lindblad master equa- 
tion, including the Kraus representation Q, the Lie al- 
gebra approach]^, the perturbative expansion 11 1 , 
the approach based on the damping bases [lj] and the 
effective Hamiltonian approach [13| . The key idea of the 
effective Hamiltonian approach is to map the Lindblad 
master equation to a Schrodinger equation by introduc- 
ing an ancilla, this leads to the advantage that almost all 
methods developed to solve(or analyze) the Schrodinger 
equation can be borrowed to solve(or analyze) the mas- 
ter equation. Leveraging on this advantage to define and 
calculate the geometric phase, as well as to formulate the 
dynamically decoherence-free subspaces [l5| is one goal of 
this paper. 

The Markov approximation is inadequate for many 
physical phenomena. Even if the environment is large 
compared to the system, it might still react on the sys- 
tem with memory. For example, the system can only 
couple to a few environmental degrees of freedom for 
short times, resulting in a memory effect of the envi- 
ronment on the system. In fact, short time scales in 
experiments often show environmental memory effects, 



a decay that can be partially undone by exploiting envi- 
ronmental memory effects [14j in the case of spin-echoes is 
an example. Also, non-Markovian quantum effects may 
play a role in the energy transfer in photosynthesis [16[. 
Hence modeling non-Markovian open quantum systems 
is crucial for understanding these experiments. There 
are many extensions developed to go beyond the Markov 
approximation]!! [H, GJ, 0, EH, [H H3. Among 
them, Breuer[24j has derived a non-Markovian master 
equation by u sing the correlated projection superopera- 
tors technique [25l [2ol [27l [2H [29| . this master equation 
can be written in a generalized Lindblad form and thus 
it is local in time. The other goal of this paper is to 
extend the effective Hamiltonian approach to this gener- 
alized Lindblad master equation, an extended damping 
basis for this non-Markovian dynamics is also given. A 
connection between these two approaches is presented. 
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This paper is organized as follows. In Sec. II, after in- 
troducing the effective Hamiltonian approach, we present 
a definition for the adiabatic and non-adiabatic geomet- 
ric phase for the open system, a connection to the ear- 
lier work is established. Based on the effective Hamil- 
tonian approach, we formulate the dynamically stable 
decoherence-free subspaces, necessary and sufficient con- 
ditions for the dynamically stable decoherence-free sub- 
space are provided and discussed. This effective Hamil- 
tonian approach is extended to a non-Markovian case 
in Sec. III. An example to calculate the effective Hamil- 
tonian and solve the non-Markovian master equation 
is given in Sec. IV. We apply this generalized effective 
Hamiltonian approach to analyze the decoherence-free 
subspaces and define the adiabatic evolution for this non- 
Markovian dynamics in Sec.V. Finally, we conclude our 
results in Sec. VI. 
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II. GEOMETRIC PHASE AND 
DECOHERENCE-FREE SUBSPACES IN THE 
EFFECTIVE HAMILTONIAN APPROACH 

The effective Hamiltonian approach [3(j is a method 
to solve the Lindblad master equation. The main idea of 
this method can be outlined as follows. By introducing 
an ancilla, which has the same dimension of Hilbert space 
as the system, we can map the system density matrix 
p(t) to a wave function of the composite system (system 
+ ancilla). A Schrodinger -like equation can be derived 
from the master equation. The solution of the master 
equation can be obtained by mapping the solution of the 
Schrodinger -like equation back to the density matrix. 
Assume the dimension of the Hilbert space for both the 
system and the ancilla is N, and let \E n (0)} and |e m (0)) 
denote the eigenstates for the system and the ancilla, re- 
spectively. The mathematical representation of the above 
idea can be formulated as follows. A wave function for the 
composite system in the iV 2 -dimensional Hilbert space 
may be constructed as 

N 

p(f)-|*(*)> = Yl Pmn(t)\E m (0))\e n (0)), (1) 
m,n— 1 

where p mn (t) = (E m (0)\p(t)\E n (0)) . Note that = 
Tr(p 2 ) < 1, i.e. this pure bipartite state is not normalized 
except when the state of the open system is pure. With 
these definitions, the master equation in Lindblad form 
(h — 1 hereafter) Q 

j t P = -i[H,p}-~Y / {LiL k p + pLiL k -2L k pLl} 

k 

= -i[H,p} + £ DP = £p, (2) 
can be rewritten in a Schrodinger -like equation [3] 

i^|*(*))=«T|*(t)>, (3) 

where TLt is the so-called effective Hamiltonian and is 
defined by 

H T = H-H A +iJ2 L £ L ki (4) 

k 

with TL = H — | ^2 k Lt^Lk- l~t A and L A are operators for 
the ancilla defined by 

(e m (0)|O A |e„(0)) = (E n (0)\O^\E m (0)). (5) 

In Eq.([2]), H is the free Hamiltonian of the system, L k are 
the Lindblad operators, and C is the Lindblad superoper- 
ator. The eigenoperators of the Lindblad superoperator 
£ compose the damping basis fl2|. namely a damping ba- 
sis vector satisfies the following eigenequation, 

LA V = \ V A V . (6) 



Note that the eigenvalues might be complex because the 
Lindblad superoperator is not Hermitian in general. One 
can also define the left eigenoperator for the Lindblad 
superoperator, 

B V C = \ v B Vl (7) 

which constructs the dual damping basis. The damping 
basis and its dual have the same eigenvalue and satisfy 
the orthonormal condition 

Tr(^B„)=<V (8) 

Now we establish the relation between the damping basis 
and the eigenstates of the effective Hamiltonian. Defining 

N 

\R») = E A» mn \E m (0))\e n (0)), (9) 

m,n=l 

where A v mn = (E m (0)\A u \E n (0)) . It is easy to check that 
Eq.© can be rewritten as, 

-M T \Ru) = \v\Rv). (10) 

This means that the pure bipartite state \R V ) defined by 
the damping basis A v through Eq.© is an eigenstate of 
the non- Hermitian operator —iTix with eigenvalue X v . 
For the dual damping basis B v , the relation is similar, 

N 

(L v \ = £ (e m (0)|(£; n (0)|^„, (11) 

m,n— 

and 

(L u \(-iH T ) = (L V \X V . (12) 

This indicates that the left vector {L v \ defined by the 
dual damping basis is the left eigenstate of the non- 
Hermitian operator —iHt- So the orthogonal and nor- 
malized condition for the damping basis can be repre- 
sented in terms of the left and right eigenvectors of the 
effective Hamiltonian, 

(L u \R^)=S uli . (13) 

To shed more light on this connection, we examine the 
steady state solution of Eq. @ in terms of both the damp- 
ing basis and the eigenstates of the effective Hamiltonian. 
The steady state means -^p = 0, which leads to Cp = 0, 
i.e. the steady state is a damping basis vector of the 
Lindblad superoperator C with zero eigenvalue. Note 
we assume C time- independent. In terms of the effec- 
tive Hamiltonian, this reads H,t\R) = 0\R). The steady 
state corresponds to the right eigenvector of the effective 
Hamiltonian TLt with zero eigenvalue. Note the existence 
of such a solution for the master equation requires that 
the determinant of the effective Hamiltonian Tlx is zero. 

With these knowledge, we now define the adiabatic 
and non-adiabatic geometric phase, and formulate the 
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dynamically decoherence-free subspaces [HI by this effec- 
tive Hamiltonian approach. To define a non-adiabatic 
geometric phase for an open system, we define a dynam- 
ical invariant operator X — X(t) for the open system by 



i^-X-[H T ,X] =0, 



(14) 



where H.t is the effective Hamiltonian of the open sys- 
tem. This definition has the same form as the dynamical 
invariant for a closed system, however the invariant op- 
erator X for open system is not Hermitian in general. 
Assume a right(left) basis {|rj}} ({|^)}) spanned by the 
right(left) eigenstates of I exists, the wavefunction of the 
composite system can be expanded in the right basis, 



l*>=£c,-(t)|» 



(15) 



By insetting Eq. (fTS"|) into the Schrodinger -like equation, 
we obtain 



(16) 



Now we prove that for non-degenerate eigenstates 
|r,->, (i = l,2,...) of I, 



(U\Ht - ij^rj) = 0, 



(17) 



for all i ^ j. As \rj) is a right eigenstate of X, we have 

l(t)\r j )=X j (t)\rj), (18) 

where Xj (t) denotes the eigenvalue corresponding to the 
eigenstate \rj). The derivative of Eg . (fl~8|) with respect to 
time yields(A J = Xj(t)), 



(k\X\rj) = XjSij + (Xj - Xi){U\fj). 



(19) 



On the other hand, from the definition of the dynamical 
invariant X, we have 

i(l m \X(t)\r n ) = (A« - X m )(l m \n T \r n ), (20) 

Eq.fdU) and Eq.® together yield, 

d 



X j 8 ij = (X i -X j ){h\[iH T + -)\r j ) (21) 



leading to, 



<M(W-i|) |r#> = 



(22) 



for i j and Aj ^ Xj. Observe that for i ^ j, we get 
Xj = 0, Vj which indicates that the dynamical invariant 
has time-independent eigenvalues, similar to the case of 
closed systems. With this result, Eq. lfll))) reduces to, 



(23) 



integrating Eq. (f23|) . we get 

Cj (t) = cj^e-'fo^^^'e-fo^^ 1 '. (24) 

This result tells us that when the system is initially in 
a right eigenstate of X, the system will remain to that 
eigenstate up to a phase factor at later times. The first 
exponential in Eq. (|24p gives the dynamical phase, while 
the second exponential generates the geometric phase. In 
a cyclic evolution, 



T (i 3 \-^h)dt'. 



(25) 



This definition was shown [3 1| to be gauge invariant and 
hence it is observable. For noncyclic evolution, a term 
arg[(Zj(0)|r,(T))] has to be added to 7*, ensuring that 
geometric is phase gauge invariant. The limitation of the 
present definition is that we have assumed the Jordan 
blocks of X are one dimensional and the eigenstates are 
nondegenerate. Beside, the success of our representation 
relies on the existence of the dynamical invariant X. For 
discussions beyond these limitations, we refer the reader 
to [13, HH. Observe that the phases defined above are 
for nonadiabatic evolution, since no adiabaticity require- 
ment has been imposed in any step of the derivation. The 
adiabatic condition can exhibit in the dynamical invari- 



ant) 



0, 



ant X(t), for a slowly varying X(t) satisfying 
we obtain \HT,X(t)] w 0. This implies that X(t) and Ht 
have a common basis of eigenstates. Thus the geometric 
phase defined above may be expressed in terms of eigen- 
states \Rj(t))(j = 1,2,...) ofT(t) as 



1 (Lj\^\Rj)dt. 



(26) 



This is the definition for geometric phases in an adiabatic 
evolution (30j. 

Decoherence remains the most important obstacle to 
experimental realizations of quantum processors. One 
well-developed method of counteracting the effects of 
decoherence is to encode quantum information into 
decoherence-free subspaces (DFSs) [M HI IM 113, H, 
l39l ]. DFSs have been defined as a set of states that un- 
dergo unitary evolution in the presence of couplings to 
the environment. However, unitary evolution of a quan- 
tum state can arise in a number of ways and this fact has 
resulted in different but related definitions for DFSs in 
the literature. In the context of Markovian master equa- 
tion, DFSs have frequently been defined as a collection 
of states p(t) fulfilling Copit) — with decoherence su- 
peroperator Co given in Eq.(f2]). Most recently, the con- 
cept of DFSs has been extended and a dynamically stable 
decoherence-free subspaces (DDFSs) has been defined [15] 
by 



d_ 

Of 



Tr[p 2 (f)] = V< > 0, with Tr[p 2 (0)] = 1, (27) 
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As shown in Ref. Eq.(f27|) is a sufficient and necessary 
condition for the quantum state p(t) to evolve unitarily. 
The definition Eq. (|2?|) together with Eq.© yield, 



d_ 

di 



Ti[p 2 (t)]=2(C D p(t)). 



(28) 



Clearly, the condition for DDFSs is less restrictive than 
that for DFSs. In the context of effective Hamiltonian 
approach, we have ^Tr[p 2 (i)] = then a 

state in DDFSs should fulfill 

d 



Of 



(*(t)|#(t)) = 0, Vi > 0, with (*(0)|*(0)) = 1,(29) 



where \^(t)) denotes the wavefunction of the composite 
system (system+ancilla) defined by Eq.([T]). We now show 
that a space spanned by 

HDDFS. = {\*l),\*2),...,\$ m )} (30) 

is a DDFSs if and only if the basis vector fulfill, 

L k \$>i) =/3 fe |$0 (31) 

for all I = 1, 2, 3, m and for all L^, as well as i?DDFSs 
is invariant under 



H 



DDFSs 



H -H J 



(32) 



k k 

By using the Schrodinger -like equation, we get from 



|<*(i)|*(t))=0 



(33) 



that 



m\ £ 4^* + 4 f 4 - £(44* + h.c.)i^)> = o. 

fe it 

(34) 

We now prove that L k \<&{) = f3 k \^i) is equivalent to 
($i|L fc ^ = (${|/9fc, where is a basis vector in 

HddfSs- By definition, any basis vector |$;) 6 Hddfss 
fulfills = /3fe| < I ) i), this equation can be rewritten as 

£|£ m (0))(£4Q)|(£ m (0)|Z fe |£„(0))|$i) =0 fc |*,).(35) 



|$;) can be expanded in the basis of the composite system 
as, 



|$,)=E*ri^(0)>|e fl (0)) ) 



(36) 



P,9 



where is the expansion coefficient. Substituting 
Eq.® into Eq.^ and noting <Ff 9 * = $f, we get 

($l\L^ = In the last derivation, we have 

used Eq.©. With L k \^i) = /3 fc |*i) and ($z|Lf = 



<$,|/3 fc , we have fTr[p 2 (t)] = (*(*)!£* 4^* + 

if 4 - _E fc (44 f + H.c.)|*(t)) = 0, where 
is an arbitrary state in TLddfSs that can be writ- 
ten as = Ei c K*)l $ i)- Thus h ddfss is a dy- 
namically decoherence-free spaces. Note that with 
Lk\^i) — the Schrodinger -like equation becomes 
*3jl*(*)) = KddfSs I *(*)>, so H DDFSs being invari- 
ant under the Hamiltonian TLddfSs is as required. To 
prove that the conditions are necessary, we suppose that 
§l{^{t)\^{t)) = holds. In general L k | can be writ- 
ten as = |*)) 



(37) 



with l^fe) denotes some state that is orthogonal to the 
state Insetting this equation into Eq. (|3"4")) yields, 



o. 



(38) 



Therefor any state |*(t)) fulfilling gj<*(*)|*(*)) = cer 
tainly satisfies Eq. ([3"Tj) . so Eq. (|3"Tj) together with Eq 
are the necessary condition for DDFSs. 



III. EXTENSION OF THE EFFECTIVE 
HAMILTONIAN APPROACH TO 
NON-MARKOVIAN DYNAMICS 

In this section, we extend the effective Hamiltonian 
approach to a non-Markovian dynamics governed by 
the generalized Lindblad master eauation|24l. [25j | . the 
corresponding generalization for the damping basis is 
also given. As mentioned above, the generalized mas- 
ter equation is obtained by the projection superoper- 
ator technique [1^, The form (Markovian or non- 

Markovian) of the equation crucially depends on the ap- 
proximation used in the derivation, reflecting in the cho- 
sen projection superoperator. When we project the state 
of the total system (system plus environment) into a ten- 
sor product, we obtain the Markovian master equation, 
while a non-Markovian master equation can be obtained 
when we use a correlated projection. In the latter case 
the master equation derived is in the generalized Lind- 
blad formdl, 



^ = -i[H k ,p k }+J2(R x kj p j R 



./A 



1 



(39) 



where H k are Hermitian operators and are arbitrary 
system operators depending on the system-environment 
interaction. If we only have a single component ps = Pi, 
this equation obviously reduces to the ordinary Marko- 
vian master equation, whereas in general cases, the state 
of the open system is ps = £ fc Pk with Trp^ < 1. 

To simplify the derivation, we first rewrite Eq. (f3"9"]) as 



i—jjT = rikPk - PkH k 



'^2 R kjPj R 



(40) 
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where Hk = Hk 



'J2j\ R jk R jk- 



Next we map this an ancilla labeled by A, 



equation into a Schrodinger -like equation by introducing 



J 



^2{E m {0)\HkPk\E n (0))\E m {0))\e n (0)) - £ (E m (0)\p k \E p (Q)){E p {0)\Hl\E n (0))\E m (0))\e n (0)) 



+*E E <i?m(o)|i?^ft|^(o))<£; P (o)|< t |i? n (o))|i : ; m (o))|e„(o)) 

n jA n 

(ft* - W£)|*k(t)> + »2^i2fcil*i(*)>, 



(41) 



where Rfcf 



and 



are operators of the auxiliary 
system defined by Eq.©, and |^ / fc(t)) is the non- 
normalized wave function corresponding to pk, defined by 
(t)|£ m (0))|e„(0)) with 

Pk,mn — 

(E m (0)\pk\E n (0)) . Eq. (|4T|) can be rewritten in a compact 
form, 



d 



i-|*(t))=H|*(t)>, 



(42) 



where |*(t)) = |%(*)), ■ • •] is a wave function 

vector, and H is an effective Hamiltonian that we will 
refer to as effective Hamiltonian matrix hereafter. The 
matrix elements of this effective Hamiltonian are, 



H 



5k j (Hk 



nt) 



(43) 



Note that the elements of the effective Hamiltonian ma- 
trix are operators and in general not non-Hermitian. The 
diagonal element includes a term Hk — H%, which de- 
scribes the free evolution for the open system and the an- 
cilla, and an interaction term i ^ A Rkk^-kk^ which results 
in the quantum jump on \^k)- The off-diagonal elements 
represent the coupling between \^k) (k — 1, 2, 3, ...), lead- 
ing the quantum jump from one component (say \^i}) to 
the other component (say |vl/j))[4(|]. When we have only 
one component, the results reduce to Ref . [l3T| . Assuming 
the Hamiltonian of the original system is time indepen- 
dent, we can obtain formally the evolution of the wave 
> 

function vector as 



!*(*)) 



iHti 



*(0)>- 



(44) 



By mapping this solution back to the density matrix 
through |* fc (t)) - E£,»=iP*,mn(t)|£m(0)}|e n (0)) with 
Pk,mn = (E m {0)\pk\E n {0)) , we obtain the time evolution 
governed by Eq. . 

Similarly, we can define a damping basis for 
the generalized Lindblad equation as a set of 



operators {A^A^,-"} satisfying —i[Hk, 



The 



definition for the dual of this damping basis is similar. 
The connection of this damping basis to the eigenstates 
of the effective Hamiltonian matrix will be given through 
an example in the Appendix. 



IV. EXAMPLE 

In this section, we will use the model and the master 
equation given in Ref.[25[ to illustrate our method. Con- 
sider a two-state system coupled to an environment. The 
environment consists of a large number of energy levels 
which constitute two energy bands. The lower energy 
band contains N\ levels while the upper one N2 levels. 
A detailed description for this model can be found in 
[ill. |43| . In the interaction picture, the non-Markovian 
master equation reads, 



(2) / 1 \ - 72 r 

Ps ( t > a _ yt cr 



ipf{t)=W-pf{t)* + 



dt 



71 r - 

2 1 



,P ( s\t)}- (45) 



where — 1,2) depend on the system-environment 

couplings as well as the energy gap of the environment. 
ct ± are the Pauli operators. Defining Hi = J2 ni \ n i)( n i\ 
and TI2 — ^2 n2 1^2) (^2 1 5 ni+n2 = Ie, where the index n\ 
labels the level of lower energy band, and n 2 is the level 
index of the upper band. The two unnormalized density 
matrixes can be obtained by pg = TrE{^iPr),i = 1,2, 
where pr is the density matrix for the whole system (the 
system plus the environment). The reduced density ma- 
trix for the system is then given by p — pg + p^ . Equa- 
tion (|43|) can be written in the form of Eq. by setting 
Hi = 0, (t = 1,2), R n = R22 = 0, #12 = ^na+, 
and i?2i = ^/72 c^ ~ ■ By the effective Hamiltonian ap- 
proach, this equation can be easily solved and an ana- 
lytical expression for the density matrix can be given as 
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follows. The elements of the effective Hamiltonian matrix where r + and t~ are the Pauli operators for the auxil- 

for Eq. are, iary system. The elements of the time evolution operator 

U corresponding to this effective Hamiltonian can be ob- 

H n = --iry 2 ( a + a - + t+t~), tained by simple algebras, 

H12 = iji o- + r+, 
H 2 i = iy 2 <7~T~ , 

H 22 = --ili{a-cT + +t~t+), (46) 



71 +72 

U 12 = 7i(l-e-^) g+T+ 

7i + 72 
72 (l_ e -(7!+72)t) _ _ 
U21 = ; cr T , 

71 + 72 

1 7 , e -(7i+72)t 

U 22 = (t + (j-t + t- + e"5 7l *(cr + cj"r"r + + a^cr+T+r") + f f cr"cr + T"r + . (47) 

7i + 72 

With this evolution operator, we can solve the Schrodinger -like equation and map the solution back to the density 
matrix, then obtain 

(1) I (7i+72e-^i+^» t )p< 1 1 ) (0)+7i(l- e ' (71+T2)t )pi 2 2 ) (0) e _i 72tp (l)( ) ^ 



PS = 



( 2) _ / Pn(0), e-i-n*($(Q) \ 

Ps = \ g-i-ntpg)^), 72(i-e-^+^' t )piV(0)+( 72 + 7ie -^+^> t )pg'(0) ] ■ ( 48 ) 



71+72 



For an initial state that only the lower band is populated, 
this solution reduces to the well-known result given in 
Ref . [25[| . We now turn to the steady state solution of 
the master equation. The existence of steady states re- 
quires that the determinant of the effective Hamiltonian 
matrix is zero. This requirement is met for H in Eq. (|46[) . 
As shown in Sec. II, the steady states are given by the 
eigenstates of the effective Hamiltonian matrix with zero 
eigenvalue. It is threefold degenerate in this example, and 
the three degenerate eigenstates in terms of the damping 
basis are 

A 01 = {p$ ) =0,P?> = \e)(e\}, 

A 02 = {pU = \g)(g\,pW = }, (49) 

A 03 = = -^\e){e\,pf = -^\ 9 ){g\) , 
I 71+72 71+72 J 

where |e) and \g) are the excited state and the ground 
state for the two-state system, respectively. This steady 
state solution tells us that there are three types of equi- 
librium for the whole system. A 01 describes an equi- 
librium that the two-state system is in its excited state 
while the environment occupies the upper band; A 02 rep- 



resents an equilibrium in which the two-state system is 
in its ground state with the environment in its lower 
band; steady state A 03 indicates that the ratio of the 
population in the ground state to the population in the 
excited state for the system is exactly the same as the 
ratio between the populations of the environment in its 
lower and upper band. One can see from these steady 
states that this non-Markovian master equation indeed 
accounts for the memory effect of the environment, lead- 
ing to an environment-state-dependent evolution of the 
system. 

V. APPLICATIONS: ADIABATIC EVOLUTION 
AND DYNAMICALLY STABLE 
DECOHERENCE-FREE SUBSPACES 

In this section, we will apply the effective Hamiltonian 
approach to define an adiabatic evolution for the non- 
Markovian open system, and present an analysis for the 
dynamically stable decoherence-free subspaces. The adi- 
abatic theorem[43] is one of the oldest and most useful 
tools in quantum mechanics. It tells us that if a state is 
an instantaneous eigenstate of a sufficiently slowly vary- 
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ing Hamiltonian H(t) at one time, then it will remain 
close to that eigenstate up to a phase factor at later times, 
while its eigenvalue evolves continuously. The notion of 
adiabaticity has been extended to Markovian open sys- 
tems in[30j, [32J, however no extension for the adiabaticity 
from Markovian dynamics to non-Markovian dynamics 
can be found in the literature. 

Having the effective Hamiltonian Eq. (|4"3")l . such an ex- 
tension is straightforward and the adiabaticity can be 
formulated as the following. An open system govern by 
the master equation Eq. VS9j) is said to undergo adiabatic 
evolution if the composite system govern by the effective 
Hamiltonian H Eq. |J3|) ( it may depends on time directly 
or indirectly) evolves adiabatically. Following the proce- 
dure in[3(|, we arrive at the adiabatic condition 



(L m (t)\R n (t)) 



An 



Ar. 



« 1, 



(50) 



where (L m (t)\ and |R n (i)) are the left and right eigen- 
state of the effective Hamiltonian matrix, respectively. 
We would like to address that the present discussion is 
restricted to systems where the effective Hamiltonian ma- 
trix is diagonalizable with nondegenerate eigenvalues. A 
generalization to the non-diagonalizable and degenerate 
case can be made by using the methods in Ref. [33| • 

To check whether the adiabatic condition Eq. (f50|) can 
guarantee the adiabatic evolution in non-Markovian open 
systems, we simulate the dynamics governed by Eq. (|45|) 
and compare it to the result from the adiabatic evolution 
defined above. To this end, we assume that the two pa- 
rameters 71 = 71 (t) and 72 = 72 (t) are time-dependent. 
In fact, the time-independent parameters 71 and 72 can 
be obtained under special conditions [29j. We define the 
following function 



(L m (t)|Rn(t)) 



(51) 



to characterize the violation of the adiabatic condition, 
where max is taken over all m and n, i.e. all the left 
and right eigenstates of the effective Hamiltonian matrix. 
The detailed expressions for the left and right eigenstates 
are given in the Appendix. We plot T as functions of 
7i(T) and ^71 (T) in Fig[T] To compare the numerical 
solution p with the adiabatic evolution p a , we use the 
fidelity[45] as a measure to quantify the difference be- 
tween the two density matrices. The fidelity is defined 
as T{p, Pa) = Tr-i/ ^fppa\ff>- It reaches one when the two 
states are the same. In Figl^l we plot the fidelity (1 — J 7 ) 
as functions of 71 (T) and -^71 (T). We choose A 03 as the 
initial state. The other parameters chosen in these figure 
are 72 (T) = 1 and J^CO = 1. Comparing FigE]with 
Fig[TJ we can find that the character function T and the 
fidelity have a very similar appearance, i.e. the smaller 
the character function is, the slighter the difference be- 
tween the numerical simulation p and the adiabatic result 
p a is. This verifies the adiabatic condition. 



1 . 



0.5- 




FIG. 1: Character function F versus 71 (T) and Jj7i (T) at a 
fixed time point T. 71 describes the decay rate of the system, 
whereas J771 characterizes the change rate of 71. Both 71 
and Jj7i can vary independently, for example when 71 = cut, 
Jj7i = uj- The other parameters chosen are 72 = 1 and 



£72 = 1. 



0.06, 




FIG. 2: The fidelity {l-T) as functions of 71 and ^71. The 
other parameters are the same as in Fig. [T] 



In the context of effective Hamiltonian approach, we 
now formulate the dynamically stable decoherence-free 
subspaces for the non-Markovian dynamics governed by 
Eq. (f3"fj]) . In terms of the wave function vector \*&(t)) — 
[|*i(t)), |*2(*)), • • - ] T , the definition of DDFSs given in 
Eq. (|2"Tj) can be expressed as 

d 



E^(^(*)l^-(*)) = o, Vi>0, 

with 5^(0)1^(0)) = 1. 

i,3 



(52) 



By the Schrodinger -like equation Eq. (|42]) , this definition 
immediately follows, 



ijk 



0. 



(53) 



If R^m) = /3fc|$m) holds for all X,k,j and the 
eigenvalue (3^ is independent of m and j, we have 
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E W ^(*i(*)l*i(*)> = 0, namely H' DDFSs = {|$ m )} 
composes a DDFSs. With R^\^ m ) = @k\®m), the el " 
ements of the effective Hamiltonian matrix reduce to 

pA pA 

+ iJ2^Rkf- (54) 

A 

Thus that Ti-' DDFSs is invariant under the effective Hamil- 
tonian with matrix elements H^P FSs is required. We 
observed that these conditions are similar to that in Eqs 
(|31|32|) . so by repeating the same procedure, we can prove 
that these conditions are necessary for DDFSs. 

VI. CONCLUSION AND DISCUSSION 

Based on the effective Hamiltonian approach, we have 
presented a self-consistent framework for the analysis of 
geometric phases and dynamically stable decoherence- 
free subspaces in open systems. Comparisons to the ear- 
lier works are made. A connection of this effective Hamil- 
tonian approach to the method based on the damping 
basis has been established. This effective Hamiltonian 
approach is then extended to a non-Markovian case with 



the generalized Lindblad master equation. As an exam- 
ple, the non-Markovian master equation describing a dis- 
sipative two-level system has been solved by this method. 
An adiabatic evolution has been defined and the corre- 
sponding adiabatic condition has been given based on 
this extended effective Hamiltonian approach. A nec- 
essary and sufficient condition for the dynamically sta- 
ble decoherence-free space is also presented. The effec- 
tive Hamiltonian approach can be extended to all master 
equations which are local in time. The geometric phase 
defined through the effective Hamiltonian is in fact a dif- 
ference of two geometric phases, when the system under 
consideration is a closed system with pure initial states. 
The present analysis is available for all master equations 
that are local in time, i.e., at any point in time the future 
evolution only depends on the present state and not on 
the history of the system. We restrict ourself in this pa- 
per to consider the case where the effective Hamiltonian 
is diagonalizable and its eigenstates are nondegenerate. 
The situation beyond this limitation can be analyzed by 
introducing Jordan blocks in the Hamiltonian, which be- 
yond the scope of this paper. 

This work is supported NUS Research Grant WBS:R- 
710-000-008-271 and NSF of China under grant Nos. 
60578014 and 10775023. 
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Appendix: Eigenstates and The Corresponding 
Eigenvalues of The Effective Hamiltonian Matrix 



and their dual 



B m = {Bf = Q,Bf = \e){e\}, 
B° 2 = {Bf = \g)(g\,Bf=0}, 
B™ = {B? = \e){e\,B? = \g){g\}, 
B^ = {Bf=0,Bf = \e)(g\}, 
B* = {B?=0 ) B? = \g)(e\}, 
B^={Bf = \e){g\,Bf = 0}, 
B 22 = {B 22 = \g)(e\,B 22 = 0}, 

72 



B 3 ={ B * = — Z!!_|e)<e|,Bi = VbXffl 



7i + 72 



7i + 72 



The alternative expressions can be given in the Hilbcrt 
space spanned by flee), \eg), \ge), \gg)} 



Af = 


= [0,0,0,0] T 


Af 


= [1,0,0,0 


T , 




= [0,0,0,1] T 


Af 


= [0,0,0,0 


T 

1 


Af = 




71 ,0,0,0 


T 

A 03 - 

! A 2 ~ 


0,0,0, 






.7i +72 






Af- 


= [0,0,0, of 


A 11 

A 2 


= [0,0,1,0 


T 


A?-- 


= [0,0,0,0] T 


A 12 

A 2 


= [0,1,0,0 


T 


Af = 


= [0,0,1,0] T 


Af 


= [o,o,o,of , 


Af = 


= [0, 1,0, 0] T 


Af 


= [0,0,0, of , 


Ai = 


[1,0,0, of, 


A\- 


= [0,0,0,-1 
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l T 



7i + 72 



In this section, we list the eigenstates and the corre- 
sponding eigenvalues for the effective Hamiltonian matrix 
H. The eight eigenvalues are, (threefold degenerate, 
corresponding to A 01 , A 02 , A 03 ) , — (twofold degener- 
ate corresponding to A 11 , A 12 ), —^i (twofold degenerate 
corresponding to A 21 , A 22 ), and —(71 +72) (correspond- 
ing to A 3 ). The corresponding right eigenstates are (in 
damping basis) 

A^ = {Af=0 7 Af = \e)(e\}, 

A 02 = {Af = \g)(g\,Af = 0}, 

72 



A 



03 



Af 



|e)(e|,^ 

71 + 72 

A u = {Af=0,A? = \g)(e\}, 
A' 2 = {Al 2 =0 7 Al 2 = \e)(g\}, 
A^ = {Af = \g){e\,Af = 0}, 
A 22 = {Af = \e){glA 22 =0}, 
A 3 = {Al = \e){e\,Al = -\g){g\} 



71 +72 



-\g)(g\ 



are the right eigenstates, and 



Bf 
Bf 
Bf = 
Bf = 
B\ 2 = 
Bf = 
B 22 = 



B\ 



[0,0,0,0], Bl 



[0,0,0,1], Bf 
303 



[1,0,0,0] 

[0,0,0,0] 

[1,0,0,0], Bf = [0,0,0,1] 

[0,0,0,0], Bf = [0,0,1,0] 

[0,0,0,0], Bf = [0,1,0,0] 

[0,0,1,0], Bf = [0,0,0,0] 

[0,1,0,0], Bf = [0,0,0,0] 
72 



7i + 72 



0,0,0 



0,0,0, 



7i 

7i + 72 . 



are the left eigenstates. 



